Triggering rogue waves in opposing currents 
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We show that rogue waves can be triggered naturally when a stable wave train enters a region of 
an opposing current flow. We demonstrate that the maximum amplitude of the rogue wave depends 
on the ratio between the current velocity, [/q, and the wave group velocity, Cg. We also reveal that 
an opposing current can force the development of rogue waves in random wave fields, resulting in 
a substantial change of the statistical properties of the surface elevation. The present results can 
be directly adopted in any field of physics in which the focusing Nonlinear Schrodinger equation 
with non constant coefficient is applicable. In particular, nonlinear optics laboratory experiments 
are natural candidates for verifying experimentally our results. 

PACS numbers: 



In the ocean, rogue waves are often observed in regions 
characterized by strong currents like the Gulf Stream, 
Agulhas Current and the Kuroshio Current [2] . Sev- 
eral ship accidents have been reported in these regions as 
being due to the impact with very large waves. One of 
these occurred in February 1986 to the SS Spray, which 
was travelling along the East coast of the USA. The ship 
was hit by a wave with a height of approximately 17 m 
(estimated by eyes from the deck of the ship), which was 
the second of a system of three consecutive large waves, 
commonly known as the three sisters. This particular 
wave system is usually observed in the nonlinear stages 
of the modulational instability process. Such instability 
was discovered in the late sixties independently by Za- 
kharov [3] and Benjamin and Feir [4 (an interesting and 
stimulating review on the subject can be found in [5 ). 
The theory is based on the linear stability analysis of a 
plane wave and predicts that a small perturbation may 
grow exponentially when eN > l/>/2, where e = k^A^ is 
the steepness of the plane wave, with /cq its wave num- 
ber and Aq its amplitude; N = uoq/IWI is the number 
of waves under the modulation, with uoq the angular fre- 
quency corresponding to the wave number and Al] 
the angular frequency of the modulation. 

The nonlinear stages of the modulational instability 
are described by the exact breather solutions of the Non- 
linear Schrodinger (NLS) equation [BHH]- Breathers are 
coherent structures that oscillate in space or time and 
have the peculiarity of changing their amplitudes as they 
propagate. They can grow up to a maximum of 3 times 
their initial amplitude and have been considered in var- 
ious fields of physics as a plausible object that describes 
the formation of rogue waves [9j [10] . Such solutions have 
also been observed in fully nonlinear simulations of the 
Euler equation, [11 . Breather solutions may also exist 
embedded in random waves [12] and can affect the prob- 
ability density function of the surface elevation and wave 



height distribution [131 HI] • We mention here that the 
concept of rogue waves is rapidly expanding to other dis- 
ciplines such as nonlinear optics and condensed matter 
(see, for example, [15H2Q]). 

In general ocean waves are characterized by a small 
value of eN . Wind seas, which are waves forced by the 
local wind field, have moderately large steepness but the 
spectral band-width only rarely allows for a large num- 
ber of the product eN that satisfies the instability crite- 
rion. On the other hand, swells, i.e. long crested waves 
that have moved out of the generating region, are char- 
acterized by a narrow spectrum (i.e. large N) both in 
angle and frequency, but they are not particularly steep. 
Hence, breathers (rogue waves) are fortunately rare ob- 
jects in the ocean. However, if swells enter into a current, 
their properties can change and, as we will show in this 
Letter, breathers solutions can be naturally triggered. 

Our analysis is based on the computation of a mod- 
ified Nonlinear Schroedinger equation, recently derived 
in [21 , that accounts for a current, U = U{x)^ assumed 
to be small with respect to the wave phase velocity Cp'. 
U/cp = 0{e). The derivation requires also that the cur- 
rent is a slowly varying function of the spatial coordinate, 
i.e., l/(A/co) = O(^), with A the entry length of the cur- 
rent (typical space scale over which the current changes) . 
The equation describes the evolution in space of the wave 
envelope A = A{x,t). In dimensional variables it takes 
the following form: 
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with Cg the group velocity, ko the wave number of the 
carrier wave and ctq = V^/cq, with g the gravity acceler- 
ation. We have found that equation ([T]) can be reduced 
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to the standard NLS equation with variable coefficients 
by applying the following transformations: 
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after which, equation ([T]) assumes the form 
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with AU{x) = U{x) — U{0); primes have been omit- 
ted for brevity. Equation ([T]) does not preserve the 
energy, which can be shown to change in space as fol- 
lows: E{x) = J \A\^dt = exp [-AU/cg] J \B\^dt. As the 
waves enter into a current, their wave height increases 
if AU < (opposing current) and decrease if AU > 
(co-propagating current). 

In the absence of a current, U=0, equation ([T]) ad- 
mit breather solutions, [6 , whose maximum amplitude 
reached during the evolution of the wave group is: 
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For sN < 1/a/2, wave groups are stable; for sN > 1/a/2, 
the amplification factor reaches the maximum value of 
three, which corresponds to the Peregrine soliton, ob- 
served experimentally recently in nonlinear optical fibers 

m- 

We will consider the evolution of a wave train initially 
in a region of zero current propagating into a stationary 
current characterized by an entry length A. We have used 
in our computation the following simple mathematical 
expression for the current: 
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Here x is the coordinate along which the waves propagate 
and Uo is the asymptotic value (large x) of the current 
that can be either positive or negative for co-propagating 
or opposing currents, respectively. For such a current 
field, we show a first example of the formation of a rogue 
wave in an opposing current in Figure [l] As initial condi- 
tions for our numerical experiments, we have considered 
a perturbed plane wave with e = 0.1 and N — 7^ which 
are typical values for an energetic swell; in terms of di- 
mensional quantities, we can imagine a wave system char- 
acterized by a period of 10 seconds (0.1 Hz). In absence 



of current, the wave train is stable and no modulational 
instability is observed. After 60 wavelengths of propaga- 
tion, the wave group enters into a current characterized 
by Uo/cg = —0.2 (opposing current) and A = lOA. At 
x/X = 60 the whole envelope grows in amplitude because 
the energy E{x) = J \A{x^t)\'^dt changes, in the pres- 
ence of a current gradient. Subsequently the envelope 
starts being strongly modulated, reaching appr oximately 
two times its local standard deviation, ^/E{x). The plot 
shows a clear example of formation of a rogue wave start- 
ing from a stable plane wave. We have performed a sys- 
tematic study on the dependence of the maximum am- 
plitude (divided by ^/E{x)) reached by the envelope as 
a function of the ratio Uo/cg which was varied from 0.1 
to 0.4 (stronger current may result in wave blocking phe- 
nomena and wave breaking). The results are shown in 
Figure |2] where the dots corresponds to our numerical re- 
sults. The plot indicates that the normalized maximum 




FIG. 1: Evolution in space and in time, normalized by the 
wave length and period respectively, of the wave envelope. 
The wave propagates for 60 wavelengths before entering into 
the opposing current of velocity Uo = — 0.25cp. The effect of 
the current is not only to increase the energy (amplitude) of 
the wave but also to trigger the formation of a rogue wave. 



amplitude increases with increasing |/7o|, evidencing that 
the originally stable plane wave has been transformed 
into a breather in the presence of a current. In the fig- 
ure we also include as a solid line the heuristic prediction 
based on the following equation: 
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Equation ([8| is a modification of the exact relation (|6|. 
The rationale behind this prediction is that the coeffi- 
cient in front of the nonlinear term in equation Q, once 
written in non-dimensional form, is the same as the one 
without current except for the exponential factor (see the 
coefficient f3{x) in ([5|). As shown in figure [2j the numeri- 
cal results are in excellent agreement with the prediction. 
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FIG. 2: The normalized maximum amplitude as a function of 
Uo/cg. The dots are the results from the numerical simula- 
tions and the solid line corresponds to equation (|8|. 
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FIG. 3: Evolution of the kurtosis as a function of the nor- 
malized distance x for different simulations corresponding to 
different values of AU/cg. For larger values of the AU/cg 
there is a clear indication of stronger deviations from Gaus- 
sian statistics. 



Once established the possibility that an opposing cur- 
rent may trigger unstable modes, it is of interest to under- 
stand whether the statistical properties of the surface ele- 
vation and in particular the occurrence of extreme events 
change as random wave trains propagate into the cur- 
rent. We concentrate our analysis on the kurtosis, z^, i.e. 
the fourth order moment of the probability density func- 
tion of the surface elevation estimated as = (t?^)/ (^^)^, 
where (...) stands for the ensemble average. In the pres- 
ence of a current, an analytical estimation of the kurtosis 
is a rather difficult task because of the nonlinearity of the 
problem. Therefore, we perform direct numerical simu- 
lations of the modified NLS equation with initial condi- 
tions characterized by the following bell-shape spectrum 
for the envelope A: 
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with AVt the standard deviation (the width of the spec- 
trum). The phases are considered randomly distributed 
in the [0,27r). Numerical simulations are computed on a 
grid of 1024 points and 600 realizations have been per- 
formed. The wave steepness V^/cqV^ was selected equal 
to 0.15 and AVL/gq = 0.2. The kurtosis is therefore esti- 
mated first as a time average and then the resulting value 
is averaged over the ensemble. The current is character- 
ized by xo = and A = lOA and different values ofUo/cg 
are considered. In Figure |3j the kurtosis is presented as 
a function of x/X. The figure clearly indicates a depen- 
dence of the kurtosis on the ratio between the current 
increment and the wave group velocity. This result is 
to some extent consistent with recent laboratory experi- 
ments, which investigated the evolution of mechanically 
generated, random wave fields over a partially opposing 
current [22]. In the absence of a current, as shown in [23 , 
the kurtosis depends on the square of the Benjamin- Feir 
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FIG. 4: f{AU/cg) = \n[{K-3)/{a7vBFI^/VS)] as a function 
of AU/cg. Dots are the results from the numerical computa- 
tion and the solid line is obtained by a linear fit. 



Index (BFI): 
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which is the ratio between the nonlinear and the linear co- 
efficients in the NLS equation properly written in nondi- 
mensional form see [12 . In the presence of the current, 
we assume the following dependence on the current ve- 
locity: 
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where a and b can be determined a posteriori from the 
numerical simulations. For a = 1 and 6 = 1, the kurto- 
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sis estimated through (11) corresponds to an analytical 



result that can be obtained by a quasi-gaussian approxi- 
mation under the hypothesis that the wave spectrum and 
the current field are slowly varying in space. The deriva- 
tion follows the one presented in [23^ performed in the ab- 
sence of current. We mention here that our aim is not to 
establish quantitatively the validity of the closure model 
but to understand how the opposing current infiuence the 
statistics of the waves. In order to verify the exponential 
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FIG. 5: The probability of exceedance for \A\/^/2E for the 
case of a linear system (dashed line), nonlinear system without 
current (solid line) and nonlinear system in the presence of an 
opposing current (dotted line). 



factor in ( 11 ), we first compute from the initial conditions 



the BFI^ then the asymptotic (large x) value of the kur- 
tosis is measured as an average value of the kurtosis over 
the last 20 wavelengths (see Figure [3|; at last we deter- 
mine the coefficient a from the simulation with A/7 = 
(the resulting value is a = 1.133). In Figure [4] we show 
the quantity f{AU/cg) = ln[(/^ - 3)/{aBFp7T/V3)] as a 
function of AU/cg. The points from the simulation lie on 
a straight line which indicates that the assumed exponen- 
tial dependence is consistent with our numerical simula- 
tions. In the figure we also show the linear fit, where the 
slope is 6 = 1.29. Moreover, we consider the exceedance 
probability (defined as p{x')dx' ^ with p{x) the prob- 
ability density function) for the envelope computed at 



=60 A for the case of AU/cg = and AU/cg 



-0.4. 



The probability of occurrence of rogue waves increases 
notably as the opposing current is stronger. Results are 
presented in Figure [5] where we have also plotted the ex- 



ceedance probability for Rayleigh distribution which is 
the one estimated for a linear process. In [24 the ap- 
pearance of rogue waves in an NLS equation with third 
order dispersion has been discussed in terms of the ratio 
between the nonlinear and linear part of the Hamilto- 
nian. In our case such ratio, proportional to the BFI, 
calculated for the initial condition is equal to 0.3; our sim- 
ulations are in the intermittent-like rogue wave regime, 
they appear and disappear erratically. 

Ocean swells are in general not very steep and wave 
packets are stable in terms of modulational instability. 
However, we have shown that breathers may be triggered 
when swells enter into a region of opposing current. This 
is an important result that should be kept in mind when 
ships navigate in the Gulf Stream or in the Agulhas Cur- 
rent or Kuroshio Gurrent in the presence of a opposing 
waves. Indeed, such currents may reach velocities up to 
1.5 meter per second and for a group velocity correspond- 
ing to waves of period equal to 10 second (a typical con- 
dition during storms), the ratio AU/cg is of the order of 
0.2, large enough to trigger a dangerous rogue wave. We 
underline that this is completely different process from 
the development of a caustic, a pure linear mechanism 
[25H27] . From a physical point of view, the mechanism of 
formation of rogue waves can be summarized as follows: 
an initial wave whose perturbation is stable in terms of 
the modulational instability may become unstable in the 
presence of a current because of the a shift of the mod- 
ulational instability band. The modulational instability 
of the wave thus leads to a triggering of the rogue wave. 
The results presented in this Letter, even though pre- 
sented in the oceanographic context, may also apply to 
the nonlinear optics physics where the nonlinearity is pro- 
vided by the type of the material. As shown in equation 
Q, the effect of the current is to change the coefficient 
of the nonlinear term in the NLE equation. New ex- 
periments characterized by materials that change their 
nonlinear properties in space could be easily performed 
and the predictions of the present work could be veri- 
fied. Acknowledgments K. Trulsen and Al Osborne 
are acknowledged for discussions. This work has been 
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